Weak indices and dislocations in general topological band structures 
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It has recently been shown that crystalline defects - dislocation lines - in three dimensional topo- 
logical insulators, can host protected one dimensional modes propagating along their length. We 
generalize this observation to the case of topological superconductors and other insulators of the 
Altland Zirnbauer classification, in d=2,3 dimensions. In general, protected dislocation modes are 
controlled by the topological indices in (d-1) dimensions. This is shown by relating this feature to 
characteristic properties of surface states of these topological phases. This observation also allows 
us to constrain these surface states properties, which is illustrated by an addition formula for (d-1) 
and (d-2) indices of a topological superconductor. 

PACS numbers: 



Topological insulators, a new phase of matter with 
protected surface states have attacted strong inter- 
est both theoretically and experimentally (for reviews, 
seei"— ). Recently we showed that dislocations in a 
three-dimensional topological insulator can host one- 
dimensional helical modes traveling along its length^. In 
addition, 3D topological insulators are characterized by 
four Z2 topological indices: one "strong index" and three 
"weak indices" The weak indices are related to the 
strong index in one lower dimension: two spatial dimen- 
sion, and we will also call them (d-1) indices. We showed 
that it is these (d-1) indices that controls the existence of 
one-dimensional helical modes hosted by dislocations^. 

3D topological insulator is just one example of many 
different types of topological band structures''^. Can 
these results be generalized to other topological band 
structures? Do crystalline defects in other topological 
band structures host topological bound states? Are there 
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(d-1) indices in general topological band structures? Is 
the existence of these bound states also governed by (d-1) 
indices? Our answers to these questions are positive. 

In this paper we discuss dislocations and lower- 
dimensional indices in general topological band struc- 
tures. If we require the topological insulator or supercon- 
ductor to have stable edge states (stable in the presence 
of impurities), there are generally 10 classes of themiSili 
according to the presence or absence of the time-reversal, 
particle-hole and chiral symmetries. One can classify 
these 10 classes of topological band structures in the 
presence of disorder (and in the absence of crystalline 
order) by, for example, studying the localization prob- 
lem on the surfaces^. The resulting classifications are 
labeled by the so-called strong topological indices in each 
dimension. We copy the table of general topological band 
structures from literature^ in the following: 
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Table: Complete list of strong indices of topological insulato rs (superconductors) and their physical realizations. 



What are the physical meaning of these strong in- 
dices? First they are topological invariants of the bulk 
band structures. They also characterize the surface(edge) 
states. For example, the edge states of the quantum hall 



insulator (class A in two dimension) is characterized by 
an integer, which is the number of one dimensional chiral 
modes traveling along the edge. And the surface states of 
a 3D topological insulator (class All in three dimension) 
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FIG. 1: (color online) Cut-glue technique in 3D(left) and 
2D(riglit). The surface(edge) states are also shown schemati- 
cally. B is the Burger's vector of the dislocation. 



is characterized by a Z2 number, which is the parity of 
the number of the surface Dirac nodes. 

On the other hand, as we will show shortly, this char- 
acterization is not complete in the presence of crystalline 
order. For example, the Dirac nodes of a 3D topological 
insulator (class All) can be at (0, 0) of the surface Bril- 
louin Zone, and it can also be at (0, tt). The difference of 
these two possible surface states are not characterized by 
the strong index. Instead it is the (d-1) indices (or the 
weak indices) that detemine the location of the Dirac 
nodes on the surfaces. In Refi, we showed that these 
weak indices also control the existence of the dislocation- 
hosted helical modes. We will generalize these results 
from class All to all classes and dimensions in the above 
Table, and discuss the new physics of general topological 
band structures in the presence of the crystalline order. 

Let us firstly provide some physical intuitions. There 
is a simple way to understand the existence of the 
dislocation-hosted helical modes in a 3D topological insu- 
lator, which is the cut-glue technique. As shown in Ref*^ 
(see FigH]), one can create a dislocation by cutting a 3D 
sample into two halves {z < and z > respectively) 
and then gluing them back nontrivially. After cutting, 
two surfaces parallel to the x-y plane are created. When 
gluing back, only the x > part of the two surfaces are 
glued back as before, and the x < part of the two 
surfaces are glued back after translating by one lattice 
spacing along y direction. This precedure creates a dis- 
location along X = z — line with Burger's vector along 
y direction. 

After cutting, the surface modes on the two surfaces 
can be described by Dirac equation. After gluing back, a 
mass term will be introduced to gap out the Dirac surface 
modes. In particular, it is clear that the gluing precedure 
actually introduces a tt Berry phase for the Dirac cone 
if it is at fcj, = TT. This means the mass term for this 
Dirac cone changes sign from a; > to a; < 0, namely 
the dislocation is nothing but a mass domain wall of the 
surface Dirac equation. It is well-known that such mass 
domain wall hosts topological pretected gapless modes- 
Application of such domain-wall-fermion bound states 



in condensed matter physics has been discussed in var- 
ious systems such as polyacetylen o^'^i and PbTe-type 
semiconductorsi^ii^. 

A careful study of these domain wall fermions in 3D 
topological insulators proves that B ■ M = 7r(mod 27r) is 
the condition to have helical mode inside the dislocation^, 
where B is the Burger's vector (real space vector) of the 
dislocation and Af is a momentum space vector composed 
of the (d-1) indices (also called weak indices in litera- 
tures). 

This cut-glue technique can easily be generalized to 
other topological band structures. We simply need to 
classify the surface states of a general topological band 
structure in the presence of crystalline order. We find 
that apart from the possible (d-1) indices, there are also 
possible (d-2) topological indices. In fact we will show 
that all the (d-1) and (d-2) indices can be naturally found 
by classiiying the surface states. And then we show that 
these surface states, after the cut-glue precedure, can 
host gapless topological bound states via, for example, 
domain wall fermion mechanism. 

Our main results is the following theorem: In the pres- 
ence of lattice order, the (d-1) indices of a d-dimension 
(d> 1, otherwise if d — 1 there is no (d-1) indices.) topo- 
logical insulator (superconductor) of a given symmetry 
class, is composed of the d strong indices, if any, of the 
(d — l)-dimension topological insulator(superconductor) 
in the same symmetry class. These (d-1) indices form a 
momentum space vector M and control the existence ( or 
the number) of the stable gapless modes hosted by a dislo- 
cation with Burger's vector B via the dot product B ■ M . 
At the same time, if d = 3, there are also three (d-2) 
indices composed of, if any, the d = 1 strong index in the 
same symmetry class. These (d-2) indices form a real 
space vector, and cannot be detected in a generic disloca- 
tion defect, but can be detect by the surface modes. We 
note that the (d-1) and (d-2) indices are presented as a 
mathematical results in a remarkable work of Kitaev^i^ 
In this work, we show these indices in a different and 
more physical fashion. The physical meaning of these 
indices, in particular the relation of these indices and 
the surface modes are explicitly shown. The well-known 
time-reversal symmetric topological insulator (class All) 
is just a special case of the above theorem. 

We prove this theorem by discussing the ten classes 
case by case. In the following we only show two cases as 
examples and the other cases are studied in Appendi:xE1 
These two examples are class D in two dimension and 
class Dili in three dimension. 

Another important result that we obtain is a general 
addition formula involving both the (d-1) and (d-2) in- 
dices: Eq.Q. This formula is useful when studying the 
interface-hosted bound states of two different topological 
band structures. 

(class T) d — 2) This is the class of p+ip superconduc- 
tor. The strong index vq € Z indicates the number of chi- 
ral majorana fermion modes at the edge. Are there other 
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FIG. 2: Deformations of edge modes of (class D d = 2) p+ip 
topological superconductor in the edge Brillouin zones. The 
particle-hole symmetry (P-H in figure) allow the motion of 
a pair of majorana edge modes from fc = to fc = tt. The 
parity of the number of Dirac nodes at fc = tt is a topological 
invariant . 



topological invariants of the edge states? Consider the 
edge of lattice vector ai . By studying the arbitrary de- 
formation of the edge modes, one finds that, by particle- 
hole symmetry, the parity of majorana fermion modes at 
fci = TT is also a topological invariant (see Figl2]). Let 
us denote it as mi (takes value of or 1) and study the 
cut-glue predure as shown in Fig IT] 

After cutting, let us assume that there is a single chiral 
majorana mode at fci = tt on each edge. At low energy 
these edge modes are described by effective theory at 
TO = 0: 



(1) 



where r's are Pauli matrices, "fix) — {'^l{x),jr{x)) 
is the two-component majorana fermion field satisfying 
7(x)* = 7(x), and jl{x),Jr{x) are the two counter- 
propagating edge modes living on the two edges. When 
we glue the two edges back, the two edges are tunneling 
to each other and a mass gap to is opened up. Note that 
this is the only non- vanishing mass term in this majorana 
theory. 

Because the majorana edge modes are both at fci = tt, 
the dislocation gives a tt Berry phase to the Dirac equa- 
tion and the mass to in Eq. ([T|) forms a domain wall at the 
dislocation. Namely 'm{x) < when x < and m{x) > 
when X > 0. Solving Eq.(IT|) one finds a single majorana 
mode 7o ~ J dx -^a "^^^ '^'^^ — jRix)) bound to 

the dislocation. One can easily check [i?, 70] — 0; i.e., 
this mode is at zero energy. Note that the majorana edge 
modes at fci = do not see the dislocation as a mass do- 
main wall and thus do not cause majorana bound state. 

Because two majorana bound states can mix and open 
up an energy gap, we conclude that oi-dislocation has a 
majorana zero mode bound to it if and only if the parity 
of majorana edge modes at fci = tt is odd. Similarly 
there is another invariant m2 at edge 122 controlling the 
existence of majorana mode in dislocation with Burger's 
vector a2. toi,TO2 are the (d-1) indices. 

In summary, the (d-1) indices can be written as a mo- 
mentum space vector M = (wiGi -|- TO2G2)/2, e Z2. 
And the condition of the existence of dislocation majo- 
rana zero mode is given hy B-M — 7r(mod27r). Although 
we only prove this condition for i? = oi , 02 , the result is 



obviously true for general B. This is because a general 
-B-dislocation can be thought as a combination of several 
dislocations with B = 0,1,0,2- 

We only studied the edge states of the a 1,02 edges. 
What about the edge states of a general edge along 
a'l = uiai -|- U2d2 direction? Another way of phrasing 
the problem is that: is M depends on the choice of the 
real space basis 0,1,0,2"! In fact we already proved that the 
answer is negative. This is because Af can be detected by 
dislocations with a general Burger's vector, and cannot 
depend on choice of basis. 

Here we comment on the relation between the current 
work and previous studies. Kou and Wen studied the 
translation symmetry protected superconductors and 
they find that there are four Z2 indices corresponding to 
the parity of the number of electrons on even by even, 
even by odd, odd by even and odd by odd lattices respec- 
tively. These four Z2 indices are related to our results and 
can be understood as follows: they are composed of the 
parity of the strong index, the two (d-1) Z2 indices, and 
one (rf — 2) =0 dimensional index. The last d = index, 
as mentioned in the footnote, is nothing but the parity 
of the total number of filled bands. This index can be 
modified by simply adding a core level, and cannot be 
detected by any low energy probe. Therefore we do not 
study this d = index in the present study. 

(class Dili d — The strong index vq ^ Z indi- 
cates the number of majorana Dirac nodes at the surface. 
Apart from vq, there are other topological invariants of 
the surface states in the presence of crystalline order. 
Consider the surface spanned by lattice vector di and 
02- By considering arbitrary deformation of the surface 
modes, one finds that, by time-reversal symmetry, the 
parity of the Dirac nodes at (7r,0), (0,7r) and (7r,7r) are 
also topological invariants. We thus have three extra in- 
variants at surface spanned by di — d2 . We re-group the 
three Z2 invariants as: the parity of the number of Dirac 
nodes along fci = tt: mi (mi = or 1), the parity of the 
number of Dirac nodes along k2 — tt: m2 (m2 = or 1), 
and a third index 7;i2. The definition of V12 is trickier: 
we choose a new surface momemtuni space basis {G'l, G2} 
such that G'l = (toiGi + TO2G2) (if mi = TO2 = then 
{G'ljGj} can be arbitrary.). We define the parity of the 
number of Dirac nodes at G2/2 as V12 (wi2 = or 1). One 
can easily see that this definition of V12 is independent 
of the choice of G2 as long as G'l , G2 form a momemtum 
space basis. We will show that TOi,m2 are (d-1) indices 
and V12 is a (d-2) index. 

By considering the cut-glue procedure in Figll] one 
finds that the parity of Dirac nodes along fc2 = tt con- 
trols the existence of majorana helical modes (a pair 
of counter-propagating majorana modes which respects 
time- reversal symmmetry.) in a dislocation with Burger's 
vector B = 02- This can be proved as follows. Assume 
there is a single majorana Dirac node at (fci, fc2) = (0, tt). 
After cutting, the low energy effective theory of the sur- 
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face modes is 

H = J dxdyj'^[{-idx)Tx + {-idy)TzfJ-z]l + m{x)^'^Tzfiy^, 

(2) 

where 7 = (xltVl, XRi Vr) ^-re the continuum hmit of the 
majorana fermions hving close to (0,7r) {L,R represent 
the left and right surfaces). Pauli matrices t is in the 
X, r] space and n is in the L, R space. Mass term m{x) 
represents the tunneling between the two surfaces. Un- 
der time reversal transformation, x ~^ ^'^'^ V ~^ ~X 
required by Kramer's theorem, namely 7 — > iTyj. Pauli 
principle requires that the matrices in the kinetic term 
is symmetric, and the matrix in the mass term is anti- 
symmetric. Finally H is time reversal symmetric which 
indicates Eq.Q is the only correct low energy theory up 
to convention. 

These Dirac majorana fermions see the dislocation as 
a TT flux tube after gluing back, because they are at mo- 
mentum A:2 = TT and the Burger's vector B = 0,2- This 
means that the mass term changes sign at the dislocation: 
mix) > when a; > and m{x) < when x < 0. Such a 
mass domain wall will have a single pair of counter prop- 
agating one-dimensional majorana fermion modes bound 
to it. This can be shown readily by solving Eq.Q as 
follows. 

Because ky is still good quantum number, we first solve 
Eq.([2]) at ky = 0. We find two zero energy bound states: 
"y^ky = 0,x) ^ -fo '>^^ 0^, where are the two 
eigenstates of Ty^y with eigenvalue +1. Because Ty^y 
and the ky matrix r^fiz commutes, and can be diagonal- 
ized simultaneously into [1,1,-1,-1] and [1,-1,1,-1] 
respectively, we immediately know that (j)± have oppo- 
site velocity along y-direction, which is the dislocation 
direction. Because these two modes are related by time- 
reversal symmetry, we name them as majorana helical 
modes. 

We only considered one Dirac node at (0,7r). In gen- 
eral, any Dirac node at fc2 = tt will contribute a pair 
of majorana helical modes. If there are even number of 
pairs, these helical modes will open up a gap by mixing. 
Thus we conclude that 7712 is odd if and only if there is 
a single pair of majorana helical modes traveling along 
dislocation with B = 0,2- Similar statement is true for 

TOl. 

We just showed that the two invariants toi, TO2 controls 
the existence of the dislocation-hosted majorana helical 
modes with Burger's vector ai, 02 respectively. Similarly 
there are four extra invariants at surface spanned by 0,2 — 
a3 and as — ai . Let us denote them by TO2 , my, (from 
02 — 03 surface) and TO3,to'j^ (from 03 — ai surface). But 
these six invariants are not completely independent. We 
must have TO2 — TO2 because they both are controlling 
the existence of majarana helical modes bound to an a2- 
dislocation. Similarly toi = m'l and TO3 = TO3. We 
thus have only three independent invariants mi, TO2, mj,. 
These are the (d-1) indices. 



We can further group them into a momentum space 
vector M = {miGi+m2G2+m3G3)/2, rrn S Z2, and the 
condition of the existence of dislocation majorana helical 
mode is B-M = 7r(mod27r). Note that we only prove this 
condition when B = oi , 02 , 03 , but the result is obviously 
true for general B. This is because a general dislocation 
can be viewed as a combination of dislocations with B = 
di, 0,2,(13. This also proves that M is independent of the 
choice of real space basis ai , 02 , 03 . 

Now we are going to establish the fact that V12 a (d-2) 
index that is related to the d = 1 strong index. Imag- 
ine we add an extra band by stacking one-dimensional 
Dili chains along a3-direction, and over ai, 02 directions. 
The edge mode of such one-dimensional Dili chain is 
a Kramer's pair of majorana modes Xi V- These edge 
modes form a two dimension lattice. Under time re- 
versal transformation: X ^ V s-i^d V ^ ~X- We can 
combine them into a complex fermion c — x + ^V- Un- 
der time reversal c — > ic^ and — > —ic. Now we can 
write down all possible lattice fermion bilinears allowed 
by time-reversal. The hopping terms are not allowed 
as tejcj t*CiC^j — ~{tc\cj)^ . Pairing terms are al- 
lowed as Ac-c] -> -A*CiCj = (Ac-c])'''. Thus these c- 
fcrmions must have a pure pairing hamiltonian on the 
di — 02 surface. In momentum space the hamiltonian 
is H = J2k ^i^)'^k'^-k + Pauli statistics requires 

that A(A:) = — A(— fc) and thus pairing must vanish and 
hamiltonian generally forms Dirac nodes at the four time- 
reversal centers: (0,0), (7r,0), (0,7r) and (7r,7r). These 
four Dirac nodes have total chirality zero and thus do 
not modify the strong index in three dimension. And 
they also do not modify mi and TO2. They only modify 
Vi2- Because the chains are along 03 directoin, they also 
do not modify TO3, V23 and V31. 

We just showed that V12 can be independently modified 
by adding an extra quasi-one-dimension band. Similarly 
V23 and V31 on the other two surfaces can also be modified 
independently. Thus there are three Z2 (d-2) indices: 
V12, V23 and V31. In the following we show these (d-2) 
indices can be summarized by a real space vector: V — 
{vi2 • 0,3 + V23 ■ di + V31 ■ a2)/2. And for a surface whose 
normal direction (a momentum space vector) given by G, 
the condition to have odd number of Dirac nodes at G2 
point on this surface is G • ^ = 7r(mod 27r). 

The above statement can again be understood by the 
stacking limit of ID Dili chains. We merely need to show: 
if we have ID Dili chains along the V = (wi2-a3+W23'ai + 
V31 ■ 02)72 direction stacked and form a 3D lattice, then 
for an arbitrary surface labeled by normal vector G, the 
chain will terminate on it with a odd number edge modes 
in a surface unit cell if and only if G • = 7r(mod 27r). 
This is obvious, as schematically shown in FigO And 
if and only if there are odd number of edge modes in a 
unit cell, when dispersion along the surface are allowed, 
they form the four Dirac nodes at (0,0), (7r,0), (0,7r) 
and (tt, tt). An even number of edge modes in a unit cell 
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two edge modes 
per unit cell 

-J 



Normal vector 
G = 2Gi + G2 



one edge modes 
per unit cell 
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Normal vector 
G = Gi + G2 



012 



a2 ^ 

2 ^ 
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FIG. 3: Stacked ID topological superconductor chains along 

V = ^ direction. If the normal vector G of a surface satisfies 

V - G = TV (mod 27r) (right) or V ■ G ^ (mod 2tv) (left), 
there will be odd or even number of edge modes per unit cell 
respectively. 



will simply open up a full energy gap when tunneling 
is allowed. These Dirac nodes explains the contribution 
to, for instance, Vi2 from this quasi-lD band. Therefore 
for a band structure with both (d-1) index M and (d-2) 
index V, one can view it as a combination of two band 
structures: a band structure with M only, and a quasi-lD 
band with (d-2) index V only. 

If one combines two general band structures together 
and ask about the total topological indices, naively one 
would guess (Ml, Vi) + {M2,V2) = (Mi + M2,Vi+V2) 
(the strong indices would obviously add up.). But this 
turns out to be incorrect. The addition of (d-2) indices 
actually involves a cross product of the (d-1) indices: 



(Mi,l/i) + (M2,V2) 
--{Mi + M2,Vi + V2^ 



, Mi X Ma 
' (2^)^ 



(3) 



where we are in the convention that the real space unit 
cell volume ai ■ 0,2 ~ 1 and the momentum space 

unit cell volume is Gi • G2 x G3 = (Stt)"^. The factor 
2 is because the (d-2) index ((d-1) index) is one half of 
a real(momentum) space lattice vector. Eq.ijS]) can be 
proved in two steps: First physically there must exist 
an addition rule of indices that makes the (M, V) space 
an additive group. And one can easily check that the 
addition rule of Eq.dH} indeed makes {M,V) space an 
additive group. This means that to prove Eq.([3]), we 
merely need to show Eq.([3]) is valid for a set of generators 
of additive group {M,V), which one can in turn check 
case by case explicitly by deforming the surface modes. 

Here we simply prove Eq.([3]) for a special case: the 
two band structures both have f = 1 , and (Mi ,Vi) — 
(^-,0), (M2,^2) = (^,0), and the combined band 
structure should have {AI,V) — ( '^i+'^a ^ m) as required 
by Eq.([3]). This is schematically shown in Fig|31 The 
readers can easily check the validity of Eq.([3]) for the 
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V12 = 1 



V23 = 
Gat 



Gi 



vr.i = 

FIG. 4: (color online) Addition of two band structures with 

{Mi,Vi) = (%-,0) and (Ala, 1^2) = (%,0) based on study 
of surface states. The shaded areas are the surface Brilloin 
Zones, and the red dots are the surface Dirac nodes. The 
resulting band structure has M = '^'-"^'^^ and V = 



more general cases. 

Eq.([3]) is general and true not only for class Dili (see 
Appendix E)) . It can be very useful when studying the 
topological bound states to the interfaces of two topo- 
logical band structures, such as a grain boundary of a 
topological superconductor. If the two topological band 
structures are labeled by {Mi, Vi) and (M2, V2), one can 
imagine to add an auxilary band structure with indices 
(-M2, -V2) to both sides, which will leave the {M2,V2) 
side topologically identical to vacuum. Such an auxilary 
band structure will not modify the interface modes topo- 
logically. We thus conclude that the interface modes is 
the same as the surface of a band structure whose indices 
are (Mi, ^i) - (M2, V2) = (Mi - M2, V1-V2- 2^^) . 

This concludes our discussion of dislocation, (d-1) and 
(d-2) indices in a general topological band structure. 
Upon completion of this manuscript, we note that a study 
by Teo and Kane^'^ of general discloations in topologi- 
cal band structures. YR thanks helpful discussion with 
Ashvin Vishwanath, Hong Yao and Dung-Hai Lee. YR 
is supported by the startup fund at Boston College. 



Appendix A: General (d-1) (and (d-2)) indices and 
dislocation-hosted bound states 



(class A d = 2): This is the well-known integer quan- 
tum hall states. The edge states are characterized by the 
number of chiral modes, which is an integer (the strong 
index). As one can easily see, there is no extra topolog- 
ical invariant as one deforms the edge states smoothly. 
Therefore there is no (d-1) index. 

(class A d = 3): This is the three dimension inte- 
ger quantum hall state. We knew that 3D quantum hall 
states can be understood by a stacking limit of two di- 
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FIG. 5: Examples of deformations of surface(edge) modes of 
topological insulators(superconductors) in the surface(edge) 
Brillouin zones, (a) All d = 3: The parities of the numbers of 
Dirac nodes at fcx = t and ky — tt are topological invariants, 
(b) All d = 2. A single Dirac node is allowed to move from 
k — to k = IT. There is no (d-1) topological invariant. 



mension quantum hall layers and are labeled by a recip- 
rocal lattice vector G that is nothing but the stacking 
direction of the layersi^. Here we pretend that we do 
not know this result and insist on classifying all possible 
surface states. We will see that we arrive at the same 
result. However, as we already showed in the main text, 
our approach is a systematical way of finding topological 
invariants of band structures and can be easily general- 
ized to other classes where new results can be obtained. 

The topologically stable surface states are the chiral 
fermi surfaces. Given a choice of the surface Brillouin 
Zone whose basis vectors are Gi and G2, the gapless 
states on this surface are characterized by two integers 
i>i and 1^2 ■ v\ is defined by the number of times that 
one crosses the chiral fermi surfaces as one moves from 
(0, — tt) to (0, tt) by tuning hi of the surface momemtum. 
If the fermi velocity of the chiral fermi surface is positive 
as it is crossed, it is counted as 1; and if the fermi velocity 
is negative, it is counted as —1. 1^2 is defined in a similar 
way by tuning k\ and count the number of times the 
chiral fermi surfaces are crossed. v\ , V2 are topological 
invariants as one cannot change them by deforming the 
surface states smoothly. And one can easily see that there 
is no more invariants, namely any two surfaces states 
with the same 1^1,1^2 can be deformed into each other. 
In addition, we can always make a new choice (denoted 
as primed) of the surface Brillouin Zone such that G'l is 
along vi ■ Gi V2 ■ G2 direction and therefore i/j = 
and v[G'i — vi ■ Gi + V2 ■ G2- These chiral fermi surfaces 
can then be understood from decoupled layered of atoms 
stacking in the G'2 direction on the surface, where each 
layer of atoms has the v'^ number of chiral edge modes. 

Now if we introduce an edge dislocation whose Burger's 
vector B parallel with this surface and terminate on 
this surface. It is clear that this dislocation must host 
n = B ■ v[G'i / {2-k) number of topologically protected chi- 
ral modes. This is because the edge dislocation can be 



simply viewed as the edge of an extra layer of atom with 
edge modes mentioned above. 

At this moment it seems that the three surfaces each 
have two integer valued topological invariants. But these 
six invariants are not independent. One can imagine 
creating a single dislocation whose Burger's vector par- 
allel with both xy-surface and yz-surface. The num- 
ber of chiral modes in this dislocation can be computed 
from either a;y-surface or yz-surface, and both results 
must be the same. This gives a constraint between this 
pair of surfaces and dictates the 1^2 on the two sur- 
faces must be identical. There are totally three such 
constraints between the three pairs of surfaces. Thus 
there are three independent topological invariants. These 
are the three (d-1) indices. And they can be sum- 
marized by a three dimension momentum space vector 
M = vi ■ Gi + 1^2 ■ G2 + 1^3 ■ G3. The number of chiral 
modes hosted by a dislocation is n — B ■ M/(27r). 

(class All d = 2) This is class of the Kane-Mele 2D 
Z2 topological insulator—. The edge modes are classi- 
fied by a single Z2 index vq which dictates the parity of 
the number of Dirac nodes. To show there is only one 
topological invariant one can simply make two ob- 
servations. First vq must be independent of the choice 
of edge. There are many ways to see this. For example 
consider a magnetic-superconductivity domain wall^i^i. 
It is easy to show that this domain wall hosts a single 
majorana fermion mode if and only ii vq ~ 1. But as the 
domain wall is moved from one edge to another smoothly, 
the majorana mode cannot disappear. Thus is inde- 
pendent of edge. Second one can see that there is no 
extra topological invariants on a given edge. For exam- 
ple, the time-reversal protected Dirac nodes must be at 
either fc = or A; = tt. One may guess that the parity 
of the number of Dirac nodes at fc = or /c = tt them- 
selves are topological stable. But one can smoothly move 
a Dirac node from fc = to A: = tt as shown in FigjS] 

(class All d = 3) This is the time-reversal symmetric 
topological insulator in three dimension and its (d-1) in- 
dices has been proven in literatures'^^''. Here we rederive 
these results using the approach of classifying the surface 
modes. 

Consider the surface spaned by lattice vector oi and 
02 of the 3D topological insulator. The parity (even or 
odd) of the total number of Dirac nodes on this surface 
is determined by the strong index. In the presence of lat- 
tice structure, the surface Brillouin zone is well defined. 
By considering arbitrary deformation of the surface fermi 
surfaces, one finds that, apart from the parity of the total 
number of Dirac nodes vq, the parity of the Dirac nodes 
along ki = IT and k2 — tt are also topological invariants 
(see Figl2) . We thus have two extra invariants at surface 
spanned by 01 — 02. Similarly there are four extra in- 
variants at surface spanned by 02 — 03 and 03 — oi . But 
these 6 invariants are not completely independent. 

To see their dependence, we consider their consequence 
in a dislocation. Consider a skew dislocation of Burger's 
vector B = 0,1. Following our discussion in literature^. 
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this skew dislocation can be obtained by the cut-glue pro- 
cedure by a cutting di — 0,2 surface containing the dis- 
location. Atoms on one side of the surface are displaced 
by _B = ai and reconnected. The Dirac nodes at fci = tt 
will see the dislocation as a mass sign-changing domain 
wall. This dislocation hosts helical mode if and only if 
the parity of Dirac nodes at fci = tt is odd. 

At the same time, we can consider the cut-glue pro- 
cedure to create same dislocation by cutting a ai — 03 
surface through it. We conclude that the parity of Dirac 
nodes at fci = tt on this ai — 03 surface must be the 
same as the one on the di — 02 surface, otherwise there 
will be inconsistency on the existence of helical mode 
inside the dislocation. Thus there is one constraint be- 
tween di — d2 surface and ai — (13 surface. Similarly 
there are another two constraints between the other two 
pairs of surfaces. Finally we conclude that there are 
three (d-1) indices which can be summarized as a vec- 
tor AI = (miGi -I- 7712(72 + 'm^Gs)/'^, iTii G Z2, and the 
condition of the existence of dislocation helical mode is 
B ■ M = 7r(mod27r). 

(class D d = 2) See the main text. 

(class D d = 3) There is no strong index and similar to 
class A, the topological protected surface modes are chi- 
ral (majorana) fermi surfaces. Given a surface spanned 
by di,a2, we can also classify the chiral majorana fermi 
surfaces by the two integers vi , V2 by the number of times 
the fermi surfaces are crossed as one tune /c2 and ki from 
— TT to TT respectively. We can also choose a new basis 
{G[,G'2} of the Brillouin Zone such that 1^2 = and 
iy[G[ = viGi + V2G2- And again the surface modes can 
be viewed as the stacking of layers of atoms along the G'l 
direction, the edge mode of each layer is v'-^ chiral majo- 
rana modes. But here, in addition to those topological 
invariants, because we are dealing with a superconductor 
with particle-hole symmetry, the parity of the number of 
chiral majorana fermi surfaces cross G'2/2 is also a topo- 
logical invariant. We denote this parity by a Z2 number 
V12 ■ One can easily see that there is no extra topological 
invariant on this surface. Similarly one can define i^'s on 
the other surfaces and V23 and Vi2. 

The ly's are the (d-1) indices and they control the num- 
ber of chiral majarana fermion modes hosted by a dislo- 
cation. Imagine an edge dislocation terminate on surface 
di — 02 with Burger's vector B = d'l. Because the sur- 
face modes can be viewed by the stacking limit, and the 
edge dislocation can be simply viewed as the edge of an 
extra disc of atoms, we immediately conclude that this 
dislocation host i^'i number of chiral majorana modes. In 
general an dislocation whose burgers vector lies within 
this surface hosts B ■ (^iGi -I- 1^2(^2) / {^ir) majorana chi- 
ral modes. 

Now imagine a dislocation whose Burger's vector B — 
di that connects the surfaces oi — 02 and di — 03. The 
numbers of chiral majorana modes can also be computed 
by the two surface must be identical. Thus the ^1 on oi — 
a2 and ai —03 must be same. One can therefore conclude 



that there are totally three (d-1) indices which can be 
summarized into a single momentum space vector M = 

■ Gi + 1^2 ■ G2 + V3 ■ G3. The number of chiral majorana 
modes hosted by a dislocation is n = B ■ M /{2tt). 

The ti's are the (d-2) indices. In the following we show 
that they are related to the d = \ strong index and form a 
real space vector V — {vi2 ■ 03 -t- ^23 • di + V31 ■ a2)/2. One 
can study an extra quasi-one-dimensional band, whose 
one-dimensional chain direction is along V. The edge 
mode of a single one-dimensional chain is a majorana 
mode X- When these chains are stacked together, the 
majorana modes can have dispersion on the surface. The 
most general lattice model of the x's, in the momen- 
tum space, will he H — e{k)xl.Xk where e(fc) is real. 
Because Xk = X^fci the Pauli statistics requires that 
e{k) = — e(— fc). Thus the surface modes must form fermi 
surfaces (zero of the real function e(fc)) that crosses the 
inversion centers (0,0), (tt, tt), (0,7r) and (tt, 0). Such 
surface mode will not have the (d-1) index, but will have 
non-trivial (d-2) index. One can also easily see that the 
(d-2) indices of a surface labeled by its normal vector G 
is non-trivial if and only if there are odd number of majo- 
rana edge modes per unit cell, which means nothing but 
G-y = 7r(mod27r). 

We thus proved that the (d-2) indices can be under- 
stood as contributed from a quasi-lD band and the (d-2) 
indices is nothing but the real space ID chain direction. 
Now the 3D class D band structure's (d-1) and (d-2) in- 
dices can be grouped as (M, V). And it can be viewed as 
the combination of two band structures (M, 0) -I- (0, V). 
It is important to note that the surface states of (M,0), 
whose band structure is the same as the stacking of 2D 
layers without ID indices along M, gives the identical 
(d-1) indices and vanishing (d-2) indices on any surface 
as {M,V). Thus the (d-2) indices is completely con- 
tributated by the (0, V) part, which is nothing but the 
stacking of ID chains. We thus proved the (d-2) index of 
an arbitrary surface is indeed given by G • V /n. 

We also find that in this case the rule of combining two 
general band structures is similar to Eq.Q: 

HMi+M2,Vi + V2 + l^^^^), (AI) 

where the only modification from Eq.Q is the factor 1/2. 
This is simply because, in our convention, the (d-1) in- 
dex M is a full reciprocal lattice vector in the present 
case. Thus adding together two (d-2)-index-absent band 
structures together may generate a non-trivial (d-2) in- 
dex. The proof of Eq. (jAl|) can also be carried out by 
studying case by case of the a few generators of {M, V) 
space. 

(class Dili d = 2) The strong index i/q ^ Z indicates 
the parity of the number of majorana Dirac nodes at the 
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edge. Consider the edge spaned by lattice vector ai. By 
considering arbitrary deformation of the edge modes, one 
finds that, by time-reversal symmetry, the parity of the 
Dirac nodes at fci = tt is also a topological invariant. And 
a dislocation of Burger's vector B — ai has a Kramer 
pair of majorana zero modes if and only if this parity is 
odd. Similarly the parity of Dirac nodes at ^2 = tt for 
the 0,2 edge is also a (d-1) index. In summary, the (d- 
1) indices can be written as a momentum space vector 
M — (miGi -|-m2G2)/2, e Z2. And the condition 
of the existence of dislocation majorana Kramer pair is 
given by _B • Af = 7r(mod27r). 

(class Dili d = 3) See the main text. 

(class AI d = 2, 3) There is no topologically protected 
edge(surface) modes. 

(class C d = 2,3) This class represents a supercon- 
ductor with full SU (2) spin rotation symmetry but no 
time-reversal. When d — 2 this is the spin quantum hall 
effecl— (basically the spin quantum number plays the role 
of charge in a usual quantum hall effect). Spin rotation 
symmetry requires even branches of the chiral mode on 
the edge. One can easily check that there is no extra 
topological invariants. So when d = 2 one only has strong 
index n where 2n is the number of branches of chiral edge 
modes. 

When d = 3 we should classify all possible chiral fermi 
surfaces that respect particle-hole symmetry and spin ro- 
tation symmetry. Similar to the discussion of 3D integer 
quantum hall, one can convince oneself that the (d-1) in- 
dices is given by three integers and can be summarized 
into a momentum space vector M = viGi+V2G2 + v^G^. 
Here, for example, 2vi is the number of times that the 
chiral fermi surfaces are crossed as the momentum is 
tuned from k2 — — tt to fc2 = tt on the ai 02 sur- 
face. The band structure can be viewed as the stacking 
of 2D spin quantum hall layers along M direction. The 
number of branches of chiral modes hosted by a dislo- 
cation with Burger's vector B is 2 ■ B ■ M, because the 
dislocation can be viewed as the edge of a disk of extra 
atoms, and the factor 2 is from spin rotation symmetry. 

(class CI d = 2,3) This class is superconductors 
with full SU{2) spin rotation symmetry and time-reversal 
symmetry. In d = 2 one can easily sees that there is no 
topologically protected edge modes and thus there is no 
topological band structures. 

In d = 3, as shown in Refi^, there are topologically 
protected surface Dirac nodes characterized by an integer 
I'Q. And the positions of these Dirac nodes are free to 
move. Thus there is no extra topological invariants. 

(class AIII, BDI, CII, d — 2,3) From now on we 
discuss the chiral classes AIII, BDI and CII. The common 



features of these classes is that the strong index in d = 1 
is z/Q G Z. This is because there is a chiral symmetry C 
that is a unitary operator, — 1, and C anticommutes 
with the Hamiltonian. The zero energy states can be 
labeled by the eigenvalues of C, ±1. If there are two zero 
modes ci,C2 with the same eigenvalue, then any bilinear 
of this two fermions C1C2, c|c2... are not allowed by chiral 
symmetry because they commute with C. Imagine we 
have a ID chain with i/Q = 1 whose left edge is labeled 
by eigenvalue -1-1 of C. Then a ladder of n-chains will 
simply have n copys of zero modes and they do not mix. 
This ladder has 1/0=11. 

The same argument holds if we stack such ID chains 
to form 2D or 3D structures. For example, if we stack 
these ID chains into a 3D crystal and have a surface 
cutting the ID chains. The edge modes of the ID chains 
remains at zero energy as one turn on tunneling between 
the chains because any bilinear is not allowed. Thus for 
non-interacting system the surface mode contributed by 
these ID chains are zero energy flat bands. If we denote 
the real space full lattice vector along which these ID 
chains are pointing along as V, the number of zero energy 
flat bands of a given surface with normal direction G is 
given by V ■ G/{2tt), which is just the number of edges 
per surface unit cell. 

In d = 2, all three classes have no strong index. The 
protected edge modes are just these zero energy flat 
bands. We conclude that the (d-1) indices of these classes 
are two integers and can be summarized as a momentum 
space vector M = viGi + V2G2- Here vi and V2 are the 
number of zero energy flat bands along di edge and 02 
edge respectively. And a dislocation with Burger's vector 
B host B-M/{2tt) number of zero modes. This is because 
an edge along B direction has B ■ M/{2tt) number of zero 
flat bands (for simplicity we assume that components of 
B does not have common divider, but one immediately 
sees that the result is general.). Thus the cut-glue pre- 
cedure immediately gives B ■ M /{2t:) zero modes bound 
to this dislocation. 

In d — 3, AIII and CII have strong indices charac- 
terized by the number of Dirac nodes on the surface. 
The positions of these Dirac nodes is arbitrary and there- 
fore have no further topological invariants. These three 
classes have no (d-1) indices. On the other hand, it is 
possible that the surface hosts zero energy flat bands 
contributed by the (d-2) index. Based on our discus- 
sion, the (d-2) index is a real space lattice vector V, and 
the number of zero energy flat bands of a given surface 
with normal direction G is given by V ■ G/ {2tt). 
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